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ABSTRACT 


The  efficient  use  of  a  search-surface  radar  or  sonar,  in 
which  one  or  more  targets  appear  on  the  screen  intermittently 
usually  demands  a  device  for  tracking  the  targets  auto¬ 
matically.  Such  a  device,  called  a  "track  while  scan  system", 
must  make  an  estimate  of  each  target's  instantaneous  position 
from  the  sampled-data  information  provided  by  the  radar. 

For  this  purpose,  an  a-B  filter  and  an  optimal  Kalman 
filter,  that  must  track  maneuvering  targets,  are  analyzed 
here  and  compared  in  terms  of  tracking  accuracy  for  tactical 
applications. 
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I .  INTRODUCTION 


Systems  such  as  search-surface  radar  and  sonar,  in  which 
the  input  data  arrives  intermittently,  frequently  require  a 
device  for  continuously  estimating  the  "present"  value  of 
the  input.  In  radar  terminology,  this  device  is  called  a 
"track  while  scan  system". 

In  a  more  general  sense,  the  term  "track  while  scan" 
system,  may  denote  any  system  which  estimates  the  "present" 
value  of  a  signal  from  the  "past"  sampled  values  of  the  signal, 
the  sampling  taking  place  at  regular  intervals. 

Vehicles,  without  maneuvering,  of  the  class  under  consider¬ 
ation  (such  as  aircraft,  ships,  and  submarines)  generally 
follow  straight  line  constant  velocity  trajectories.  If  the 
vehicles  were  not  able  to  deviate  from  these  trajectories, 
i.e.,  could  not  maneuver,  then  the  tracking  problem  could  be 
solved  quickly  and  simply  using  standard  filtering  such  as 
an  a-g  filter. 

Historically,  a-g  filters  were  designed  to  minimize  the 
mean  square  error  in  filtered  position  and  velocity  under  the 
assumption  of  small  velocity  changes  between  data  samples. 

Thus,  most  a-g  filters  have  little  capacity  to  track  targets 
that  either  accelerate  or  maneuver  (changing  direction) .  An 
important  early  paper  [1]  defines  6  in  terms  of  a  with  a  a 
design  parameter.  This  greatly  simplified  the  use  of  maneuver 
detectors  by  reducing  the  design  problem  to  a  single  variable. 
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A  later  advance  in  ot-S  filter  design  was  achieved  by 
minimizing  the  mean  square  error  in  predicted  position 
[2]  . 

A  system  that  developed  using  the  above  technique  is  the 

r 

VEGA  LN  I  (Thomson-CSF,  France). 

The  final  choice  of  parameters  for  the  a- 8  tracker  will 
always  end  in  a  compromise  between  smoothing  the  input  noise 
(measurement  errors)  and  retaining  some  ability  to  follow  a 
maneuvering  target.  Various  performance  measures  have  been 
used  for  this  compromise.  For  example/  steady  state  noise 
variance  reduction  and  the  ability  to  follow  a  target  capable 
of  responding  to  impulse  accelerations  are  used  as  performance 
measures  in  [1]  to  derive  good  steady  state  filter  parameters. 
Time  varying  noise  variance  reduction  and  the  ability  to  follow 
a  randomly  maneuvering  target  are  performance  criteria  which 
follow  from  [3] .  In  this  paper,  such  a  filter  is  analyzed, 
the  frequency  response,  stability,  noise  characteristics  and 
transient  error  are  derived  and  plotted,  and  a  scheme  for 
optimizing  the  two  dynamic  parameters  (>*-$)  is  suggested. 

A  criterion  for  a- 6  filter  design  with  a  numerical  and 
graphical  example  is  also  given. 

Since  the  majority  of  tactical  weapons  systems  requires 
that  manned  maneuverable  vehicles,  such  as  aircraft,  ships, 
and  submarines,  be  tracked  accurately,  an  optimal  Kalman  filter 
has  been  derived  for  this  purpose,  based  on  the  early  work  of 
[4,5].  The  target  model  for  tracking  applications  must  be 
sufficiently  simple  to  permit  ready  implementation  in  weapons 


systems  for  which  computation  time  is  at  a  premium  yet 
sufficiently  sophisticated  to  provide  satisfactory  tracking 
accuracy. 

The  target  acceleration  selected  for  the  Kalman  filter, 
and  hence  the  target  maneuver  is  correlated  in  time;  namely, 
if  a  target  is  accelerating  at  time  t,  it  is  likely  to  be 
accelerating  at  time  t+t  for  sufficiently  small  t.  By  sim¬ 
plifying  the  maneuvering  model  used  in  the  Kalman  filter  above, 
the  state  vector  can  be  reduced  from  six  to  four  elements. 

The  model  simplification  -  simplified  Kalman  filter  -  is 
achieved  by  assuming  (incorrectly)  that  the  vehicle's  change 
in  velocity  is  uncorrelated  between  samples;  i.e.,  the  maneuver 
is  white.  In  the  comparison  which  follows,  the  Kalman  filter 
and  Simplified  Kalman  filter  are  compared  with  the  a-B  filter 
in  a  variety  of  tactical  environments  with  tracking  sensors, 
utilizing  Monte  Carlo  simulation  techniques  on  realistic 
target  trajectories  to  verify  their  theoretical  performance. 


II.  FUNDAMENTALS  OF  TRACK-WHILE-SCAN 


Any  system  which  performs  a  tracking  function  must  obtain 
and  utilize  the  basic  target  parameters  of  position  and  rates 
of  motion.  In  the  earlier  system  both  position  data  and 
velocity  data  were  used,  to  maintain  the  tracking  antenna  on 
the  target  at  all  times,  thus  limiting  the  system  to  the  classi¬ 
cal  one  target-at-a-time  tracking  function. 

In  a  track-while-scan  system,  target  position  must  be 
extracted  and  velocities  calculated  for  many  targets  without 
holding  the  radar  antenna  fixed  on  one  target.  Obviously  in  a 
system  of  this  type,  target  data  is  not  continuously  available 
for  each  target  track  at  a  rate  dependent  upon  the  scan  rate 
of  the  system. 

In  a  typical  TWS  system  the  data  rate  is  one  unit  of  data 
per  second  or  a  scan  rate  for  the  search  radar  antenna. 

Since  the  antenna  is  continuing  to  scan,  some  means  of 
storing  and  analyzing  target  data  from  one  update  to  the  next 
and  beyond  is  necessary.  The  digital  computer  with  its  memory 
and  computational  capability  is  employed  to  perform  this 
function  and  also: 

(i)  To  provide  a  tactical  picture 
A  display  of  the  tracks  of  all  vehicles  observed  by  sensors 
showing  present  positions,  courses  and  speeds  etc.,  is  essen¬ 
tial  for  the  deployment  of  a  ship  and  its  weapons.  The  computer 
enables  rapid  use  to  be  made  of  sensor  information  thereby 
ensuring  that  the  picture  is  accurate  and  up-to-date. 


(ii)  For  use  in  processes  such  as  threat  evaluation 
and  weapon  assignment 

The  computer  can  forecast  future  likely  positions  of  tracked 
vehicles  and  rapidly  perform  necessary  calculations  to  assist 
operators  in  the  assessment  of  threats  and  in  the  optimum 
weapon  deployment  to  deal  with  them. 

(iii)  To  provide  target  information  fox  weapon 
deployment 

The  computer  can  be  used  to  generate  smooth  tracks  from  noisy 
information  thereby  being  able  to  pass  information  to  a  weapon 
sensor  more  accurately  than  information  derived  from  a  single 
plot  on  display.  Hence,  the  central  concept  underlying  any 
TWS  system  is  that  the  sensor  itself  continues  to  perform  its 
primary  function  of  search  (scanning)  and  data  input  while  the 
remainder  of  the  system  performs  the  target  tracking  function. 
The  sensor  function  simply  provides  target  position  data  to 
the  computer  subsystem  where  target  velocities  and  position 
prediction  are  calculated.  In  a  military  application  the 
major  advantage  of  a  TWS  system  is  the  elimination  of  the 
process  of  target  designation  from  a  search  radar  to  a  fire 
control  radar. 

The  tracking  information,  developed  in  the  TWS  system, 
is  used  as  a  direct  data  input  to  the  computation  of  a  fire 
control  solution. 

Therefore,  as  soon  as  a  target  is  detected  a  fire  control 
solution  is  available  without  the  inherent  delay  caused  by 
the  designation  process.  The  time  required  from  first  detection 
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to  fire  control  solution  is  on  the  order  of  miili-seconds  for 
a  TWS  system  as  opposed  to  seconds  or  even  minutes  for  a 
manually  designated  system  employing  separate  search  and  fire 
control  sensors. 

The  focus  of  the  following  chapter  has  been  to  answer  the 
question:  "What  functions  should  the  TWS  system  perform  in 

order  to  combine  the  search  and  tracking  tasks  into  one 
integrated  unit?" 

A.  A  TRACK  WHILE  SCAN  METHOD 

The  method  of  solving  the  track  while  scan  problem  is 
based  upon  the  assumption  that  the  radar  furnishes  target 
position  information  once  each  scan.  The  scheme  can  be  imple¬ 
mented  by  a  combination  of  special  radar  circuits  (hardware) 
and  software.  Existing  systems  also  provide  for  operator 
modification  of  system  tracking. 

Any  track-while-scan  system  must  provide  for  each  of  the 
following  functions: 

(i)  Target  detection 

(ii)  Generation  of  target  acquisition  and  tracking 
"windows"  or  "gates" 

(iii)  Target  track  initiation  (assignemnt  of  targets 
to  track  files) 

(iv)  Target  data  input  and  track  correlation 

(v)  Track  "window"  prediction,  smoothing  and 
positioning 
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1.  Target  Detection 

Target  detection,  localization  and  designation  are 
accomplished  by  the  radar  sab-system  in  the  usual  sense. 

Hence,  the  TWS  receives  the  following  data  from  the  outside, 
via  wire  links: 

(i)  The  surveillance  antenna  azimuth 

(ii)  The  radar  synchronization  (presyne) 

(iii)  The  surveillance  radar  video 

2.  Generation  of  Target  Acouisition  and  Tracking  "Windows" 


A  window"  cin  be  defined  as  a  small  volume  of  space 
initially  centered  on  a  target,  which  will  be  monitored  on 
each  scan  for  the  presence  of  target  information.  Each  initial 
target  detection  will  cause  a  relatively  large  acquisition 
"window"  to  be  generated,  centered  on  the  position  of  the 
detection.  When  a  target  is  initially  detected  the  algorithm 
receives  only  the  position  data  for  that  initial,  instantaneous 
target  position.  The  acquisition  window  is  then  generated, 
for  example: 

Initial  position  =  (Range,  Bearing,  Elevation) 

Range  window  =  R  ±  1000  yars  (914.1  meters) 

Bearing  window  =  B  ±  5°  (0.0873  radians) 

Elevation  window  =  E  ±  5°  (0.0873  radians) 


The  acquisition  window.  Fig.  la,  is  large  in  order  to  allow 
for  target  displace.,  ent  during  the  following  scans  of  the 
radar.  If  on  the  next  radar  scan  the  target  is  within  the 
acquisition  window,  a  tracking  window  is  generated  in  the 
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same  manner  as  the  acquisition  window.  Although  Fig.  lb 
shows  only  the  very  small,  120  yards,  1.5°,  1.5°  tracking 
window,  in  actual  practice  intermediate  window  sizes  are 
generated  until  a  smooth  track  is  achieved. 

3 .  Track  Initiation 

Concurrent  with  the  generation  of  the  acquisition 
window  a  track  file  is  generated  in  order  to  store  the  posi¬ 
tion  and  window  data  for  each  track.  In  addition  to  the  basic 
position  and  window  data,  calculated  target  velocities  and 
accelerations  are  also  stored  in  each  track  file.  Track  files 
are  stored  within  the  digital  computer  (or  processor)  sub¬ 
systems  memory  and  the  data  is  used  to  perform  the  various 
calculations  necessary  to  maintain  the  track. 

Each  track  file  occupies  a  discrete  position  of  the 
digital  computer's  (or  processor's)  high  speed  memory.  As 
data  are  needed  for  computation  or  new  data  are  to  be  stored, 
the  portion  of  memory  which  is  allocated  for  the  required  data 
will  be  accessed  by  the  system  'programs  (software) .  In  this 
manner  a  diversity  of  data  in  addition  to  the  tracking  data 
may  be  stored  in  the  "track"  file,  for  example,  ESM  data,  IFF 
information.  The  generation  of  the  track  file  begins  with  the 
initial  storage  of  position  data  along  with  a  code  to  indicate 
that  an  acquisition  window  has  been  established.  If  target 
position  data  is  obtained  on  subsequent  scans  of  the  radar,  the 
file  is  updated  with  the  coordinates;  the  velocities  and 
accelerations  are  computed  and  stored,  and  the  acquisition 
window  code  is  canceled. 


The  acquisition  window  is  then  decreased  in  size  to 
the  tracking  window  and  the  track  code  is  stored  which  indi¬ 
cates  an  active  track  file. 

As  the  radar  continues  to  scan,  each  input  of  data 
is  compared  with  the  window  positions  of  active  track  files 
until  the  proper  file  is  found  and  updated.  However,  it  should 
be  noted  that  the  search  for  the  proper  track  file  is  generally 
not  a  sequential  one-to-one  comparison.  This  method  is  much 
too  slow  to  be  used  in  a  system  where  speed  of  operation  is 
one  of  the  primary  goals. 

This  idea  of  comparing  output  data  with  window  posi¬ 
tions  leads  us  to  the  problem  of  correlation  of  data  input  to 
track  files  and  what  to  do  if  correlation  is  ambiguous. 

4.  Resolution  of  Track  Ambiguity 

Track  ambiguity  arises  when  either  multiple  targets 
appear  within  a  single  track  window  or  two  or  more  windows 
overlap  on  a  single  target. 

This  occurrence  can  cause  the  system  to  generate 
erroneous  tracking  data  and  ultimately  lose  the  ability  to 
maintain  a  meaningful  track.  If  the  system  is  designed  so 
that  an  operator  initiates  the  track  and  monitors  its  progress, 
the  solution  is  simply  for  the  operator  to  cancel  the  erroneous 
track  and  initiate  a  new  one. 

For  systems  which  are  to  be  automatic,  software  (pro¬ 
gramming)  decision  rules  must  be  established  which  will  enable 
the  tracking  program  to  maintain  accurate  track  files.  Decision 
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rules  as  to  target  input  data  and  track  correlation  can  be 
logically  developed  with  an  understanding  of  the  definition 
of  the  track  ambiguity  problem. 

5.  Track  Window  Prediction,  Smoothing,  and  Positioning 
In  a  track-while-scan  system  tracking  errors  also 
exist  due  to  target  motion.  The  tracking  window  now  has 
replaced  the  "tracking  antenna"  and  this  window  must  be  posi¬ 
tioned  dynamically  on  the  target  in  a  similar  manner  as  was 
the  "tracking  antenna".  However,  there  is  no  "servo"  system 
to  reposition  and  smooth  the  tracking  window's  motion.  This 
repositioning  and  smoothing  must  be  dene  mathematically  within 
the  TWS  algorithm.  To  this  end,  smoothing  and  prediction  equa¬ 
tions  (Eqs.  2. 1-2. 3)  are  employed  to  calculate  the  changing 
position  of  the  tracking  window.  Instead  of  the  system  "lagging" 
the  target  the  tracking  window  is  made  to  "lead"  the  target 
and  smoothing  is  accomplished  by  comparing  predicted  parameters 
with  observed  (measured)  parameters  and  making  adjustments 
based  upon  the  errors  derived  from  this  comparison. 

Figs.  2  a-b  and  3,  illustrate  the  TWS  principle, 
track  positioning  and  the  TWS  general  processing  loop  of  the 
alpha-beta  type,  respectively. 

B.  DEFINITION  OF  TRACKING  EQUATIONS  FOR  a- 6  FILTER 

Tracking  in  a  track-while-scan  (TWS)  system  consists 
basically  of  two  functions: 

a.  Smoothing 

b.  Correlation 
f 


22 


NOIiLISOd 


a.  OWS  PRINCIPLE 


ACTUAL  TRACKING 


FIG.  2.  TWS  PRINCIPLE  AND  POSITION  TRACKING 


FIG.  3.  TWS  GENERAL  PROCESSIN1 


Smoothing  is  the  processing  of  the  sensor  reports  to  derive 
an  estimate  of  both  target  velocity  and  position. 

Correlation  is  the  sorting  of  sensor  reports  into  groups, 
to  determine  which  belong  to  which  target  then  predicting  a 
new  coordinate  on  which  to  center  the  correlation  region. 

The  constants  of  proportionality,  a  and  6,  used  in 
correcting  the  position  and  velocity,  respectively,  of  the 
estimated  target  course,  completely  characterize  the  perform¬ 
ance  of  the  TWS  system.  These  constants  are  the  dynamic 
parameters  or  so-called  "smoothing  constants"  of  the  system. 

Finally,  the  simplest  case  target  tracks  are  based  on 
smoothing  and  prediction  of  an  alpha  (a)  -  beta  (B)  tracker 
operating  in  a  cartesian  coordinate  reference  frame.  The 
information  is  ordinarily  obtained  from  a  coordinate  converter 
operating  on  the  raw  polar  to  cartesian  transform.  The  ct-3 
filter  described  by  the  well-known  a-B  tracking  equations  [6,7] 

Xg  =  Xp  +  a (X^  -  Xp)  smoothing  Equation  (2.1) 

N  »N-1  ® 

Vg  =  Vg  + - - prediction  Equation  (2.2) 

Xp  =  Xg  1  +  ^  T  prediction  Equation  (2.3) 

where : 

Xg  =  smoothed  position 


predicted  position 


smoothed  velocity 


=  measured  position 


=  filter  parameters  or  filter  gains 


=  sampling  time  or  time  between  detections 


Substituting  Eq.  (2.3)  into  Egs.  (2.1)  and  (2.2)  yields 


Xg  =  (l-oOXg-1  +  (1-a)  TV^_1  +  axJJ 


vS  - 


-  |  Xs_1  +  +  I  XM 


Xp+1  =  x”  +  T 


(1-a)  (1-a) T 


(1-8) 


V 

L  s  J 


a  [XM)J 


=  fl  T]  X£ 


V  ! 
S  i 

J 
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These  equations  represented  in  block  diagram  form  appear  in 
Fig.  4. 


Fig.  4.  Block  Diagram  of  an  ot-6  Tracker 


where  the  following  matrices  and  vectors  are  defined: 


(1-cO 

(l-ot)T 

a 

A  = 

B  = 

-e/T 

(1-6) 

I 

6/T 

C  *  [1  T]  X  = 


1 .  Stability  Analysis 

Since  the  system  characteristic  equations  are  in  the 
form  of  difference  equations,  Z-transform  theory  is  helpful 
in  determining  filter  transfer  functions  and  stability  criteria 
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For  the  system  of  the  form: 


Vl  -  “«  +  B<XK  '  CV 

or 

Vi  -  <A  -  Bc)x„  +  bxm 


thus, 


G  (Z) 


X(Z) 

VZ) 


[ZI  -  (A-BC) ]-1  B 


(2.9) 


Using  Eq.  (2.9),  the  Z-transform  for  each  of  the  several 


components  of  Eg.  (2.7)  and  (2.8)  has  been  computed  and  is 


given  below.  Hence,  the  a-g  tracker  transfer  functions: 


Gxs{z) 


XS(Z) 

x^zT 


Z [a (Z-l)  +  8 ] 
D  ( Z ) 


Smoothed  (2.10) 
position 


GVS(Z) 


VS(Z)  _  fz(Z-l) 
x^TzT  =  D(Z) 


Smoothed  (2.11) 
velocity 


GXP(Z) 


XP(Z)  a (Z-l)  +  BZ 

XM(Z)  D(Z) 


Predicted  (2.12) 
position 


Vp(Z)  _  !<z-i) 

xM(z)  =  dTzT 


Predicted  (2.13) 
velocity 


28 


where  the  characteristic  equation  of  the  system  is 


D(Z)  =  | ZI  -  A 1  =  0 


or 


D(Z)  =  Z  -  (2  -  cx  -  8)  Z  +  (1  -  a)  =  0 


(2.14) 


The  stability  of  the  system  can  be  determined  from  the  loca¬ 
tion  of  the  roots  of  the  characteristic  equation  above. 

A  few  methods  are  available  for  determining  whether 
or  not  a  polynomial  in  Z  contains  a  root  or  roots  on  or  out¬ 
side  the  unit  circle.  One  method  is  to  modify  the  Routh- 
Hurwitz  stability  criterion.  The  Routh  stability  criterion 
tells  whether  or  not  any  of  the  roots  of  a  polynomial  lie  in 
the  right  half  of  the  complex  plane. 

Since  the  following  transformation 


Z  = 


r  +  1 
r  -  1 


maps  the  interior  of  the  unit  circle  in  the  Z  plane  to  the 
lef-half  r  plane,  with  this  transformation,  permits  the 
application  of  the  Routh  stability  criterion  to  the  polynomial 
in  r  as  in  continous-time  systems. 

Now  (2.14)  becomes  a  polynomial  in  r,  in  the 


form 


0 


(2.15) 


2 

Sr  +  2ar  +  4  -  2a  -  6  = 

So,  the  resulting  requirements  for  stability  are 

a  >  0  ,  S  >  0  and  (2a  +  3)  <  4  (2.16) 

An  additional  stable  condition  exists  when  s  =  0.  Thus  the 
resulting  necessary  and  sufficient  conditions  for  the  stability 
of  the  track-while-scan  system  are 

a  >  0  ,  6^0,  and  (2a  +  8)^4  (2.17) 

These  inequalities  determine  a  "stability  triangle"  in  the 
a-S  plane,  for  which  all  internal  points  and  all  points  on 
the  base  (8=0)  in  the  interval  0  <  a  <  2  correspond  to  a 
stable  system.  This  triangle  is  shown  in  Fig.  5. 

The  conditions  for  underdamped,  critically  damped,  and 
overdamped  transient  response  are  found  by  inspecting  the 
sign  of  the  discriminant  of  Eq.  (2.14).  The  resulting  condi¬ 
tions  are: 


(a+g)2  <  4B  </ . .  ->  underdamped  (2.18) 


8  <  1,  (a+8)2  =  46  C-  ■  >  critically  damped  (2.19) 


2 

a  <  1,  6  <  i,  (a+B)  >  46  <=>  overdaraped 


i 

i 


(2.20) 


FIG.  5.  THE  STABILITY-  TRIANGLE 


All  other  values  of 
(a ,3)  inside  the  stability 
triangle  and  on  the  base 
($  =  6)  in  the  interval 

0  <  a  <  2 


The  transient  response 
contains  at  least  one 
demped  oscillatory  natural 
mode  with  a  rate  of 
oscillation  equal  to  one 
half  the  sampling  frequency 


The  permissible  values  of  a  and  3  are  shown  in  Fig.  6 


Fig.  6.  Allowable  Values  of  a  and  0  in  a-B  Tracker 

Further  restrictions  will  be  placed  on  this  region,  when  the 
frequency  response  characteristics  of  the  system  are  examined. 

2 .  Frequency  Response  Characteristics 

In  this  section,  the  frequency  response  characteristics 
of  the  a -6  filter  are  formed.  This  approach  has  not  been 
done  before,  in  detail. 

The  frequency  response  of  the  filter  can  be  found  by 
placing  Z  =  e-|luT  and  e-|ajT  =  cos  uT  +  j  sin  o>T  into  Eqs.  (2.10)~ 
(2.13) .  These  equations  are  now  in  the  form: 


[a  (cos2ioT-oosu)T)  tScosu/T]  +  j  [a  (sin2u)T-sina/T)  +8siimT] 

D(e^) 


(2.21) 
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li 


.  m  [^(cxds2u)T-c30SuT)  ]  +  j  [|-(sin2wT-sinuT)  ] 

^  1  '  - - ^ - 


(2.22) 


<W>T> 


_  [(g+B)cosa/r  -  g]  +  j  [  (g+B)  sintoT] 
D(e>T) 


(2.23) 


[^•(cdscoT-I)  ]  +  j[|sinaiT] 


(2.24) 


where 


D(e^U)T)  =  [cos2uiT- (2-g-B)  cosaiT+ (1-a)  ]  +  j  [sin2uT-  (2-g-B)  sinuiT] 


The  amplitude  and  phase  characteristics  of  Gxg  (smoothed 
position)  and  Gxp  (predicted  position)  are  plotted  and  shown 
in  Figs.  7-16,  for  several  values  of  g  and  B.  Also  the  ampli¬ 
tude  and  phase  characteristics  of  Gvs  (smoothed  velocity)  and 
Gw  (predicted  velocity)  are  shown  in  Figs.  17-26.  All  the 
amplitudes  and  phases  are  plotted  as  a  function  of  g,  3  and 
0)T  . 

Observing  these  figures  and  Eqs.  (2.7)  -  (2.8)  one 
finds  that  Xg  is  the  result  of  passing  X^  through  a  low  pass 
filter,  is  the  result  of  differentiating  X^,  and  g  should 
never  be  larger  than  one. 

It  may  be  seen  also,  for  positional  smoothing,  if 
a  =  0,  all  sensor  information  is  ignored  whereas  if  g  =  1, 
there  is  no  smoothing  of  positional  information.  Similarly 
B  =  0  causes  sensor  information  to  be  ignored  in  the  estima¬ 
tion  of  velocity  whereas  8  >  1  will  cause  overcorr ection . 
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radians)  AMPLITUDE 


1.00  l.S 


3  =  0.5 

3  =  0.75 


FIG.  11.  AMPLITUDE  AND  PHASE  OF  SMOOTHED  POSITION  G 
FOR  a  =  0.75 


FIG.  13.  AMPLITUDE  AND  PHASE  OF  SMOOTHED  POSITION  C, 
FOR  a  *  0.9 


AMPLITUDE  AND  PHASE  OF  PREDICTED  VELOCITY  G 
FOR  a  *  0.1 
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AMPLITUDE  AND  PHASE  OF  SMOOTHED  VELOCITY  G 
FOR  a  =  0.9 
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FIG.  24.  AMPLITUDE  AND  PHASE  OF  PREDICTED 
VELOCITY  FOR  a  =  0.9 


Considering  the  figures  for  several  values  of  a  <  1, 
a  seems  to  control  the  bandwidth  of  the  low-pass  filter  and 
8  has  more  control  over  the  damping.  In  fact  8  should  be 
somewhat  smaller  than  a  such  that  resonant  spikes  do  not 
occur. 

Better  accuracy  occurs  between  smoothed  and  predicted 
positions,  when  a  has  very  large  values  compared  with  8. 

Table  I,  which  follows,  summarizes  the  accuracy  and 
provides  an  important  design  tool.  In  general  the  frequency 
of  the  input  signal,  the  sampling  time  T,  and  the  filter 
parameters  j,  and  p  control  the  filter's  response. 
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ACCURACY  (0-1)  BETWEEN  SMOOTHED  AND  PREDICTED 
POSITION  FOR  a)T  =  0  to  wT  =  TT 


X 

0.1 

0.5 

0.75 

0.9 

0.01 

l 

1.0-0.90 

1.0-0.98 

1.0-0.99 

l 

) 

0.05 

1.0-0.60 

1.0-0.90 

1.0-0.94 

0.10 

1.0-0.33 

1.0-0.82 

1.0-0.89 

0.25 

1 

1.0-0.71 

0.50 

! 

1.0-0.50 

0.75 

1.0-0.33 

ACCURACY  (0-1)  BETWEEN  SMOOTHED  AND  PREDICTED 
VELOCITY  FOR  0)T  =  0  to  (uT  =  it 


\a 

&X 

0.1 

0.5 

0.75 

0.9 

0.01 

1.0* 

1.0* 

1.0 

0.05 

1.0 

1.0 

1.0 

0.10 

1.0 

1.0 

1.0 

0.25 

1.0 

0.50 

1.0 

0.75 

1.0 

*As  we  expect  100%  accuracy,  since  the  vehicle  was 
generally  assumed  to  follow  straight  line  constant 
velocity,  when  an  a-8  filter  is  used. 


TABLE  I.  ACCURACY  BETWEEN  SMOOTHED  AND  PREDICTED, 
POSITION  AND  VELOCITY 
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III.  NOISE  CHARACTERISTICS  OF  ct-8  FILTER-ERRORS- CRITERIA 


In  the  study  of  any  filter  it  is  essential  to  know  the 
characteristics  of  the  desired  signals  and  the  noise  which 
excites  the  filter.  It  is  desirable  also  to  know  how  the 
choice  of  a  and  8  affects  the  degree  to  which  the  noise  is 
smoothed  or  exaggerated  by  the  system.  The  description  of  the 
noise  processes,  prediction  errors  and  methods  (criteria) 
for  designing  an  a-S  filter,  proceeds  as  follows: 

A.  NORMALIZED  NOISE  POWER  OF  PREDICTED  POSITION 
By  making  the  following  change  of  variable: 


AX  =  Vs  T 


(3.1) 


Equations  (2.7)  and  (2.8)  describing  the  filter  can  be 
rewritten  as: 


r- 

N 

- 

r 

N-l 

xs 

(1-a) 

(1-a)  J 

xs 

a 

i 

l 

+ 

AX 

-8 

(1-6)  : 

AX 

6 

l  - 

L 

j 

L  J 

[Xp]N+1 

=  tl  1] 

V 

AX  I 


'V 


N 


(3.2) 


(3.3) 


where  now 
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(1-a)  (1-ct) 


(1-3) 


=  (1  1] 


!  Xs  ■ 


x  = 


The  mean  and  covariance  equations  are  defined  as 


(X]N 


.'X^-1  +  B'X? 


(3.4) 


=  A,PN-1A,T  +  B*a2B'T 


(3.5) 


where 


[X]N  =  EtX11]  ,  the  expected  value. 


P  p  1 

11  *12  I 


P  = 


P  D 

21  *  22  j 


Pn  =  E^xs-xs)(xs-xs)] 


P12  =  E[(XS-XS)  (AX  -AX)]  =  P21 


P22  =  E[  (AX  -  AX)  (AX  -  AX)  ] 


ov„  =  standard  deviation  of  measurement  error. 
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The  covariance  equation  for  the  filter  becomes 


N 

>u; 

r  2 

|  (1-a; 

| 

2 (1-a) 2 

(1-a)2 

’Pli' 

1— 1 
to 

1! 

1  "6  (1-a) 

» 

i 

(1-28)  (1-a) 

(1-6) (1-a) 

P12 

■  ! 

1  i 

T*» 

i 

2 

n  t  1  o  \ 

/  1  n  \  2 

L  22, 


+  a8 


-28 (1-3) 


'4' 


(1-6) 


,  N-l 


(3.6) 


The  variance  of  the  predicted  position  is  easily  computed 
in  terms  of  the  variances  from  Eq.  (3.6) 


P11  +  2  P12  +  P22 


(3.7) 


The  steady  state  solution  of  Eq.  (3.6)  is  computed  by  leting 
P (N)  equal  P(N-l)  and  the  resulting  algebraic  equations  are 
solved; 


23  ~  3aB  +  2a' 
a  (4  -  2a  -  8) 


(3.8) 


6 (2a  -  8) 
a(4  -  2a  -  8) 


(3.9) 
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(3.10) 
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POSITION 


where : 


N 

v 


0 


2 

v 


N  M 

v  v 


additive  measurement  noise 
02  for  N  =  M 

0  for  N  ?  M 


and  by  letting 


Variance  (X^"1’’'") 


N+l 
VAR (Xp  ; 


Variance  Reduction  Factor  =  VRF 


i 


hence 


VRF (Xp+1) 


VAR (Xp+1) 


(3.14) 


or 


VRF  (X 


N+l. 
P  ' 


VAR(Xp+1) 


2g2+2  g+a  g 
a (4-2ci-g ) 


(Steady  (3.15) 
state) 


The  VRF  is  plotted  in  Fig.  28  as  a  function  of  a  and  g.  Com¬ 
paring  this  figure  with  Fig.  27  (normalized  noise  power  of 
predicted  position) ,  the  two  figures  look  alike,  but  this  is 
not  true,  when  both  graphs  are  plotted  together  in  Fig.  29. 


C.  PREDICTION  ERROR  DUE  TO  CONSTANT  ACCELERATION 

Constant  acceleration  X(t)  results  in  fixed  error  in 
predicted  position  X^+^  on  reaching  steady  state,  i.e.,  for 
N  -*■  °°  . 
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NORMALIZED  PREDICTION  ERROR  DUE  TO  ADDITIVE  NOISE 
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Specifically  for  noiseless  constant  acceleration  trajec¬ 


tories, 


(3.16) 


the  steady  state  prediction  error  becomes 


* 


e 


lim 

oo 


••  2 
XT 

8 


(3.17) 


Common  names  for  steady  state  prediction  error  e  due  to 


constant  acceleration  are: 

(i) 

Dynamic  error 

(ii) 

Truncation  error 

(iii) 

Systematic  error 

(iv) 

Bias  error 

Typical  design  procedure  for  balanced  dynamic  and  random 
errors  is 


...Standard  deviation, 
of  random  error  ’ 


.Dynamic, 
'  error 


that  is 


} 


VAR (Xp+1) 


—  e 


(3.18) 


D.  TRANSIENT  ERROR 

The  transient  error  of  the  system  becomes  significant  when 
the  input  is  switched  from  one  time-shared  target  to  another 
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or  when  the  target  makes  a  sharp  maneuver.  Since  a- 8  trackers 
aie  initialized  by  original  estimates  of  position  and  velocity 
and  in  many  applications  these  estimates  are  so  inaccurate  that 
the  filter  exhibits  transient  errors  before  settling  to  a 
steady-state  solution.  In  the  transient  phase,  when  the  esti¬ 
mates  are  poor,  one  would  like  to  rely  on  the  measurements 
and  thus  weight  the  errors  more  heavily.  After  several  measure¬ 
ments  have  been  made,  the  accuracy  in  prediction  increases. 

So,  the  transient  error  yields 


| 


i 

} 


| 

l 

c 

l 


_  r  ,  N+l  N  2  _  T2(2  -  a) 

VN+1  '  P  "  ciS  (4  -  2a -B) 

n-  u 


(3.19) 


The  transient  error  D^+1  is  plotted  and  shown  in  Figs.  30-32, 
as  a  function  of  a,  8  and  T.  Observing  these  figures,  one 
finds,  that  the  sampling  time  T  greatly  affects  the  transient 
error.  For  large  values  of  T  the  transient  error  increases. 
An  interesting  phenomenon  occurs  for  small  values  of  a  and 
small  values  of  8.  The  transient  error  increases  sharply. 


E.  OPTIMAL  STEADY  STATE  RELATIONS  FOR  THE  a -8  TRACKER 
The  simplest,  possible  system  is  one  in  which  a  and  8 
are  fixed  constants.  The  theory  of  such  systems  has  been 
studied  in  [8].  For  such  systems,  T.R.  Benedict  and  G.W. 
Bordner  [1]  applied  the  calculus  of  variations  to  derive  the 
following  relationship: 
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0.75 
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FIG.  32.  TRANSIENT  ERROR  FOR  T  =  0.5  AS  A  FUNCTION 
OF  a  AND  3 
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s 


(3.20) 


2 
a 

( 2-a ) 

A  few  years  later  S.R.  Neal  [7]  used  the  results  of  linear 
estimation  theory  to  derive  the  relation: 

(2a  +  g) 2  =  80  (3.21) 


As  a  matter  of  interest,  the  above  Eqs.  are  plotted  and 
compared  in  Fig.  33.  Observing  these  figures  it  is  inter¬ 
esting  to  note  that  the  relations  (3.20)  and  (3.21)  are  quite 
similar  for  0  <_  a  <  0.4  and  0  <_  3  £  0.1. 

Constant  parameter  systems  suffer  from  the  incompatible 
demands  that  good  smoothing  requires  heavy  damping  (i.e.,  small 
values  of  a  and  0  -  small  noise  response.  Figs.  27-28),  while 
good  resporse  to  maneuvers  requires  light  damping  (i.e.,  large 
values  of  a  and  6  -  small  transient  error.  Figs.  30-32).  Light 
damping  and  therefore  large  noise  response,  can  lead  to  low 
probabilities  of  weapon  sensor  acquisition  and  problems  of 
plot-to-track  association.  Heavy  damping,  implying  poor 
maneuver  response  (large  transient  error) ,  can  cause  sudden 
loss  of  cracks  (sometimes  termed  track  death)  through  failure 
to  associate  with  subsequent  plots. 

These  limitations  led  some  workers  to  optimize  for  varia¬ 
ble  parameter  systems,  where  a  and  0  are  varied  according  to 
the  state  of  the  track.  Some  systems  have  been  developed 
wherein  a  and  0  were  initially  selected  arbitrarily  and  changed 
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during  program  development  by  trial  and  error.  Various 
operational  sets  of  values  being  derived  for  various  states 
of  track.  Such  methods  are  adaptive  and  are  usually  economi¬ 
cal  in  computer  use,  both  in  terms  of  required  storage  and 
run  rime,  but  generally  have  no  theoretically  optimum  adaptation 

F.  DISCOUNTED  LEAST-SQUARES  CRITERION 
To  find  linear  trajectory 


f  i 


Xp  =  Xp  +  N  T  Vp 


(3.22) 


.  .  .  N 

wr.ich  minimizes  sum  of  weighted  errors  between  Xp  and 

yN  „N-1  N-2  . 

M'  ^4  '  i*e.. 


'N 


l  <XM~r  XP)2  ^ 
r=0 


(3.23) 


where  0  _<  0  <  1,  still  too  many  degrees  of  freedom  for 
selecting  gain  terms,  a  and  8.  Minimize 


{eN} 


(xJJ-Xp)2e°  +  ( x^j- X-Xp- 1 ) 2 0 1  +  (x”"2-Xp  2)02  +  ... 

(3.24) 


by  solving  the  above  equation  in  order  to  weight  the  differ¬ 
ence  between  measured  and  predicted  values.  This  yields 
in  simple  gain  terms 


a  =  (1  -  0*) 


(3.25) 


s  i 


j  i 
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6 


(3.26) 


(1  -  0) 


2 


Common  names  for  discounted  least-squares  a-6  filter  are 
(i)  Critically  damped  a-S  filter 
(ii)  Fading-memory  polynomial  filter  of  degree  1 
Equations  (3. 25) -  (3.26)  are  plotted  as  a  function  of  (theta) 
in  Fig.  34. 

More  recently,  processes  have  been  developed  in  which 
and  are  made  to  change  with  time  in  order  to  continually 
compute  the  least-squares  line  through  the  observations.  Such 
approaches  assumed  that  errors  are  equally  distributed  in 
x  and  y  and  had  a  constant  standard  deviation.  The  formulae 
for  changing  a  and  0  in  this  manner  were  worked  out  in  [9] . 

For  the  incorporation  of  the  nth  measurement: 


2  (2N  -  1) 
N(N  +  1) 


(3.27) 


8 


_ 6_ _ 

N  (N  +  1) 


(3.28) 


The  above  equations  are  plotted  as  a  function  of  the  number 
of  measurements  N,  in  Fig  35. 

It  is  clear  that,  for  large  N,  a  and  8  tend  to  zero, 
i.e.,  observations  will  be  increasingly  ignored.  This  sug¬ 
gests  that  there  should  be  some  maximum  value  of  N.  The 
maximum  value  used  is  generally  7  to  15.  Such  a  method  may 
be  made  adaptive  if  a  means  of  detecting  changes  in  motion 
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is  used,  i.e.,  if  turn  detection  is  provided.  Then,  if  a 
turn  is  detected,  the  values  of  a  anc  B  may  be  raised  simply 
by  lowering  N.  Doing  this  will  impro\  the  turn  following 
capability. 

1.  Numerical  Example 

Assume  the  following; 


Range  measurement  =  a 


30.5  ft 


max 


A 


Desire  /VAR(Xp+1) 


5g  =  150  ft/sec 

31.6  ft 


Find  suitable  critically  damped  a-B  filter  design. 
By  substituting  Eqs.  (3.25)  and  (3.26) 

a  =  1  -  e2 

8  =  (1  -  9)2 

in  Eq.  (3.15) 


VRF  = 


VAR(X^+1) 


2a  +  ag  +  23 
a (4  -  2a  -  3 ) 


and  letting 


VAR<x£+1)  =  (31. 6)2,  a2  =  (30. 5)2, 
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raa 


yields  • 


VRF  =  1.074 

0  =  0.5 


hence 


a  =  0.75 

0  =  0.25 


From  Eq.  (3.18) 


e*  =  3  J VAR (Xp+1)  =  3  J(3 1.6) 2  =  94.8  ft 

but  from  Eq.  (3.17) 

*•  2 

*  XT 

E  =  -  — 5 — 


hence 


T  =  0.3848  sec. 


Finally,  from  Eq.  (3.19) 


D  N+i 
XP 


T2 (2  ~  a) 
aB(4  -  2a  -  8) 
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yields 


D  XTJ,  =  D  =  0.439  ft2 

VN+1 

XP 

The  values  for  a  =  0.75  and  6  =  0.25  satisfy  the  critically 
damped  Eq.  (2.19) 

( a  +  3 )  2  =  43  and  8  £  1 

The  above  values  for  a- 3#  VRF  and  D  can  also  be  determined 
by  using  Figs.  27-31-34. 
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IV.  ESTIMATING  OPTIMAL  TRACKING  FILTER 
FOR  MANEUVERING  TARGET 


The  work  which  follows  is  an  attempt  to  apply  optimum 
filter  theory  to  the  tactical  radar  environment.  There  it  is 
desired  to  obtain  optimal  estimates  of  target  positions  and 
their  predicted  tracks. 

There  are  several  possibilities  for  structuring  a  digital 
filter.  From  the  work  of  Kalman  [3],  the  filter  should  con¬ 
tain  the  same  dynamic  model  as  that  of  the  incoming  signal 
shown  in  Fig.  36. 


Measurement  Noise  VK 


Excitation 

Signal  Process 

Signal 

Measurement 

Observable 

Dynamics 

*K 

Process 

yk 

i 

» _ i 

H 

a.  Signal  Process 


b.  Filter 

Fig.  36.  Schematic  of  Filter  Structure  in  Relation  to  Signal  Process 


78 


The  modeling  of  the  various  system  components  involved 
in  a  tactical  weapon  system,  such  as  tha  radar  measurement 
process  and  the  target  itself,  is  essential  in  the  design 
of  practical  tracking  and  control  algorithm.  By  modeling 
the  target  to  be  tracked  and  the  accuracy  of  the  radar's 
measurements  [10],  then  a  practical  tracking  procedure,  con¬ 
sistent  with  the  computer  limitations  and  weapon  system 
requirements,  can  be  designed. 

A.  SENSOR  AND  VEHICLE  MODELING  -  DYNAMIC  EQUATIONS 

The  tracking  systems  under  consideration  utilize  sensors 
that  provide  measurements  of  range  and  bearing.  The  selection 
is  intended  to  reflect  that  this  pair  of  measurements  is  most 
common,  however,  other  output  measurements  as  elevation  and 
range  rate  (Doppler)  are  often  available. 

The  vehicles  to  be  tracked  can  be  modeled  by  the  state 
equations : 

X  (K+l)  =  <j>X(K)  +  Gu(K)  (4.1) 

where 

Y(K) 

Y(K) 

0  (K) 

0  (K) 

=  vehicle  state  vector  at  time  t  =  KT 
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u(K)  = 


ux(K) 


u2(K) 


change  in  vehicle  range  rate  between 
time  K  and  time  K+l 

change  in  vehicle  bearing  rate  between 
time  K  and  time  K+l 


=  state  transition  matrix 


(4.2) 


G  = 


T  =  sampling  period. 


The  tracking  sensor  measures  target  position,  range  and  bearing 
or  elevation  and  provides  the  following  output  equation: 


Y  (K)  =  HX (K)  +  V  (K)  * 


where 


measured  rate  at  time  K 


measured  bearing  at  time  K 


(4.3) 


H  = 


and  V (K)  = 


VX(K) 

v2(k) 


(4.4) 
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The  measurement  noise  covariance  matrix  R(K),  satisfies 


R (K)  =  E[V(K)V  <K)3  = 


R(K) 


(4.5) 


0  (K) 


assuming  the  noise  V^(K)  and  V2OO  are  independent.  The 
selection  of  sensor  coordinates  (R,  6),  at  this  point,  for 
the  covariance  matrix  R(K)  has  been  made  because  the  output 
matrix  H  assumes  the  extremely  simple  form  shown,  and  the  R(K) 
becomes  diagonal.  When  Doppler  measurements  are  available, 
this  selection  of  sensor  coordinates  becomes  extremely  advan¬ 
tageous  because  the  cartesian  forms  for  H  and  R(K)  becomes 
complex  and  time  varying  and  often  impose  computational 
penalties  for  real-time  implementation. 

Consider  the  following  polar  to  cartesian  transformation 
(i.e.,  from  R,e  to  X,Y) 


X  =  R  sin  9 


Y  =  R  cos  6 


(4.5a) 


For  the  Kalman  filter  in  XY  coordinates,  the  measurement 
covariance  matrix  R(K)  is  a  function  of  the  radar-target 


geometry.  Therefore,  the  elements  of  the  covariance  matrix 
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2 

°XX 

R(K)  = 

2 

°XY 

are 

2  2  2  2  2  2 

ovv  o_  cos  6  +  R  aQ  sin  6 

aA  o 

2  2  2  2  2  2 

avv.  -  cr_,  sin  0  +  R  oa  cos  e  (4.5c) 

y  x  k  y 

2  2  2  2 
aXY  =  ^aR  ~  R  a0^  S^n  9  COS  9 

2  2 

where  o_  and  a.  are  the  variances  of  the  range  (R)  and 

K  u 

bearing  (e)  measurement  errors,  respectively.  It  is  inter¬ 
esting  to  note  that  in  general,  the  coordinates  after  the 
transformation  (X  and  Y)  are  not  independent.  The  singular 

cases  where  they  are  independent  occur  for  0  =  0°,  90°, 

oT 

10°,  ...  or  for  R  =  .  These  special  cases  are  easily 

explained  from  the  assumption  that  the  measurement  errors  are 
Gaussian,  with  typical  contours  of  equal  probability,  which 
are  ellipses. 

Hence,  if  the  axes  lie  in  the  directions  of  the  axes  of 
the  frame  of  reference,  the  ccvar:  ance  terms  RU’)  R^21 
are  2ero  and  the  independence  of  X,Y  errors  for  0=0°,  90°, 
utc. 

When  R  =  a_/a  _ ,  these  ellipses  reduce  to  circles  which 

K  n 

may  be  considered  as  limiting  ellipses  whose  major  and  minor 
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axes  lie  along  the  axes  of  the  reference  frame.  Finally , 
since  the  observation  y(K)  consists  of  x  and  y  position 
measurements,  hence,  the  observation  V(K)  can  be  assumed  to 
have  covariance 


R(K) 


(4 . 5d) 


B.  STATISTICAL  DESCRIPTION  OF  TARGET  MANEUVER 

The  input  sequence  u(K)  is  additive  ' .ace  (or  maneuver) 
noise  that  results  in  the  vehicle  deviating  from  a  constant 
velocity  trajectory. 

Although  the  maneuver  history  and  observation  noice  are 

T 

not  independent,  the  covariance  E[u(K)V  (j)]  is  zero.  Indeed, 
the  radar  cross  section  of  a  piloted  vehicle  changes  during 
a  maneuver,  causing  the  radrr  observation  noise  to  depend  on 
the  particular  maneuver  being  exercised  by  the  vehicle. 

The  maneuver  noise  is  neither  Gaussian  nor  white.  For 
example,  the  pilot  of  an  aircraft  moving  at  constant  velocity 
will  generally  not  maneuver  unless  threatened  by  either  radar 
detection  or  attackinc  vehicles.  His  maneuver  will  then  often 
be  a  turn  or  an  increase  or  decrease  in  his  forward  velocity. 

A  typical  maneuver  probability  density  is  shown  in 
Fig.  37. 
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Fig.  37.  Typical  Probability  Density  of  Target  Maneuver 

The  quantity  A  is  the  maximum  acceleration  which  the  plane- 
pilot  combination  can  withstand.  Values  of  the  density 
between  non-maneuver  (u  =  0)  and  maximum  (u  =  ±A)  are  non 
zero  because: 

(i)  The  vehicle  may  not  be  accelerating  a '  its 
maximum  rate 

(ii)  The  projection  of  a  circular  maneuver  on  any 
dimension  can  give  values  of  u  from  -A  to  A. 
Clearly,  then,  the  maneuver  density  is  not  Gaussian. 

C.  DISCRETE  TIME  EQUATIONS  OF  MOTION 

It  is  often  desirable  to  whiten  the  maneuver  noise  so 
that  system  equations  to  which  optimal  filtering  theory 
applies  can  be  obtained.  This  is  done,  as  in  [4]  The 
whitening  procedure  for  a  discrete  signal  is  analogous  to 
the  procedure  developed  by  Wiener  and  Kolmogorov  to  white 
continuous  signals  [11].  So,  u(k)  may  be  expressed  recursively 
in  terms  of  the  white  noise  sequence  W(K)  by 
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u (K+l ) 


pu  (K)  +  W  (K) 


(4.6) 


Using  the  above  equation,  the  following  set  of  system  equations 
are  obtained  having  white  noise  sequences  W^(K)  and  W2 (K) 
as  their  only  inputs: 


X  (K+l)  =  <f>X(K)  +  GW  (K) 


(4.7) 


Y  (K)  =  HX  (K)  +  V  (K) 


(4.8) 


where 


X(K) 


r(K) 
r(K) 
(K)  i 


u 


l 


1 

0  (K) 
0  (K) 


vu(K) 

J 


(4.9) 


and 


W(K) 


W1(K) 

w2(k) 


(4.10) 
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Also 


4>  = 


1  T  0  0 
0  110 
0  0  p  o 
0  0  0  1 
0  0  0  0 
o  0  o  o 


0 

0 

0 

T 

T 

X 

0 


0 

0 

0  j 
0 
1 

p . 


■ 

State  Transition 
Matrix 


(4.11) 


G  = 


(4.12; 


1  0  0  0  0  0 

0  0  0  1  o  0- 


[Observation  matrix] 


(4.13) 


and  where 


Q(K)  =  E[W(K)WT(K)  ] 


°  °M2(1-02> 


Covariance  of 
the  measured 
noise 


(4.14) 
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2 

The  values  for  the  maneuver  variances  a...  and  the  correlation 

Mi 

coefficient,  p,  depend  on  the  maneuver  characteristics  of 
the  vehicles  being  tracked. 

2 

Hence,  the  variances,  of  target  acceleration  are 

calculated  using  the  model  illustrated  in  Fig.  37.  The 
target  can  accelerate  at  a  maximum  rate  A  (-A)  and  will  do 
each  with  a  probability  P^  and  a  probability  P2  of  not 
accelerating  at  all,  with  an  assumed  uniform  probability 
distribution  of  amplitude: 

1  -  (2P,  +  P.) 

P(u)  =  - - - —  (4.15) 


ef  accelerating  between  -A  and  +A. 

The  acceleration  variable,  therefore,  has  mean  zero  and 
variance 


A2(l  +  4PX  -  P2) 
3 


(4.16) 


Consequently,  the  variables  u-^  and  u2,  which  are  assumed 
independent,  have  zero  means  and  u^  has  variance 


2 

°M1 


A2T2(1  +  4P1  -  P2) 
3 


(4.17) 


while  u2  has  variance 


2 

°M2 


A2T2(1  +  4PX  -  P2) 


3R 


(4.18) 
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where  R  is  the  target  range  from  the  sensor  and  where  all 
these  quantities  have  appropriate  units. 

The  correlation  coefficient  (|p|  <  1)  can  be  modelled 


E [u (K) u (K-l) ] 
P  -  2 

°Mi 


1  -  pT 


T  <  — 


T  >  — 
P 


(4.19) 


The  quantity,  p,  is  essentially  the  inverse  of  the  average 
maneuver  duration.  This  correlation  model  is  analogous  to 
that  in  [4],  which  has  the  discrete  time  form  r(K)  =  a...  e  , 

Mi 

where  a  is  the  inverse  of  the  continuous  maneuver  time  con- 
stant  of  the  target.  Hence,  p  equals  e  .  When  aT  is  small, 
p  can  be  approximated  closely  by  1-aT,  so  that  a  and  p  become 
identical . 

The  two  extreme  cases  occur  when  p  is  unity  and  when  p 
approaches  zero.  So,  the  first  case  represents  the  completely 
correlated  case,  and  the  second  is  the  completely  uncorrelated 
case,  namely  white  noise. 


D.  FILTER  DESCRIPTION 

Three  types  of  filters  are  considered  as  potential  can¬ 
didate  algorithms  for  tracking  vehicles  that  are  described 
by  the  model  just  discussed.  These  filters  are  the  Kalman, 
the  simplified  Kalman  and  an  a-8  filter. 
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1.  The  Kalman  Filter 


The  Kalman  filter  uses  the  augmented  version  of  the 
model  presented  earlier  in  order  to  obtain  white  excitation 
(maneuver)  noise. 

The  method  of  computing  the  optimum  estimate  (the  filter) 
is  as  follows: 


X(K|K) 


X(K|K-1)  =  <f-X  (K-l  j  K-l)  (4.20) 

X(K|K-1)  +  P(K|K-1)HT[HP(K|K-1)HT+R(K) i"1 


.  [y (K) -HX (K i K-l) ] 


(4.21) 


where 


P(K|K-1)  =  <j)P  (K-l !  K-l)  ij^+GQ  (K-l)  GT 


(4.22) 


P(K|K)  =  P(K|K-1)-P(K|K-1)HT[HP(K|K-1)HT+R(K)]  1HP ( K | K- 1 ) 

(4.23) 


In  these  equations 

(i)  The  double  argument  always  denotes  an  estimate 

(ii)  X(K|K)  =  Filtered  estiamte  of  X(K) 

(iii)  X (K  j  K-l)  =  one-step  predicted  estimate  of  X(K) 

(iv)  P(K|K)  =  Covariance  of  filtered  estimate 

(v)  P(K|K-1)  =  Covariance  of  predicted  estimate 

(vi)  (y (K) -HX (K | K-l) ]  =  Error  in  the  predicted 

observations 


(vii)  G(K) 


=  P(K|K-1)HT[HP(K|K-1)HT+R(K) ] 


=  Kalman  Gain  Matrix 


=  Is  a  matrix  of  adjustment  coefficients. 
The  matrix,  G(K),  reflects  the  relative 
confidence  one  should  have  in  the 
observed  data  as  compared  to  the 
predicted  estimate. 


The  filter  is  initialized  on  the  basis  of  two  observations 
as  follows : 


y1(2) 

|ty1(2)-y1(l)3 

0 

y2<2) 

|ty2(2)-y2(1)] 


(4.24) 


The  elements  of  the  corresponding  covariance  matrix,  P(2|2), 


are: 


(4.25) 

2  2  ‘ 

The  Quantities  a...,  a„.  should  be  calculated  as  discussed 

Ml  M2 

previously  (Eq.  4.17-4.18). 

A  block  diagram  of  the  discrete  Kalman  filter  is  shown  in 
Fig.  38. 


^jK  =  ^SclK-l  +  | K-l^ 

^(K-l  =  **K-l|K-l 


Fig.  38.  Block  Diagram  and  Equate  ,iis  of  the  Discrete  Kalman  Filter 


2 .  The  Simplified  Kalman  Filter 


By  simplifying  the  maneuver  model  used  in  the  Kalman 
filter,  the  state  vector  can  be  reduced  from  six  to  four 
elements,  and  the  number  of  independent  components  of  the 
covariance  matrix  from  ten  to  six.  The  model  simplification 
is  achieved  by  assuming  (incorrectly)  that  the  vehicle's 
change  in  velocity  is  uncorrelated  between  samples,  i.e., 
the  maneuver  is  white. 

The  regular  Kalman  filter  requires  two  augmented  state 
variables  in  order  to  whiten  the  target  maneuver.  If  the 
maneuver  is  assumed  white,  no  augmentation  need  be  performed 
and  the  simplification  just  discussed  occurs.  This  simpli¬ 
fied  Xalman  filter  can  therefore  also  be  referred  tc,  as  an 
unaugmented  Kalman  filter. 

The  utility  of  this  filter  is  greatest,  therefore, 
when  either  sensitivity  of  tracking  performance  to  assumed 
maneuver  correlation  is  small,  or  when  the  target  maneuver 
approaches  whiteness  relative  to  the  sensor  data  rate. 

The  equations  for  this  filter  and  all  quantities 
except  the  following  have  previously  been  defined: 


QUO 


Covariance  of 
"assumed  white" 
maneuver  noise 


(4.26) 


and  the  elements  of  the  corresponding  covariance  matrix, 
P(2|2),  are: 
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3.  Recursive  Algorithm 

The  problem  is  solved  recursively  by  first  assuming 
the  problem  is  solved  at  time  K-l.  Specifically  it  is  assumed 
that  the  best  estimate  X(K-ljK-l)  at  time  K-l  and  its  error 
covariance  matrix  P (K-l | K-l)  are  known. 

(i)  Calculate  the  one-step  prediction 

X  (K  |  K-l)  =  4>X  (K-l  |  K-l)  (4.28 

(ii)  Calculate  the  covariance  matrix  for  the  one-step 
prediction 

P(K|K-1)  =  <t>P(K-l|K-l)  <j)T  +  GQ(K-1)GT  (4.29) 

(iii)  Calculate  the  prediction  observation 


y(KjK-l) 


HX (K | K-l) 


(4.30) 


(iv)  Calculate  the  filter  gain 


G (K)  =  P(K|k-1)HT[HP (KjK-i)HT+R(K) ]_1  (4.31) 

(v)  Calculate  the  new  smoothed  estimate 

X(K|K)  =  X(K|K-1)  +  G(K) [Y(K)-Y(K|K-1) ]  (4.32) 


P(K|K)  =  [I  -  G (K) H]  P(K|K-1)  (4.33) 

In  summary,  starting  with  an  estimate  X(K-l|K-l)  and  its 
covariance  matrix  P(K-l|K-l)  after  receiving  a  new  observation 
Y(K)  and  calculating  the  six  quantities  in  the  recursive 
algorithm,  a  new  estimate  X(K|K)  and  its  covariance  matrix 
P (K | K)  are  obtained. 


V.  IMPLEMENTATION  AND  SIMULATION  RESULTS 

In  order  to  evaluate  the  three  filter  algorithms  in  a 

variety  of  tactical  environments#  an  air  vehicle  type,  and  two 

tracking  sensors  were  selected  for  analysis. 

Maneuver  statistics  (A,  P^,  P^#  u)  were  selected  to  match 

2  2 

the  vehicle,  and  sensor  statistics  (a-,,  a a ,  T)  were  selected 

k  y 

for  each  combination  of  sensor  and  data  entry  evaluated. 

One  trajectory  was  constructed  for  the  vehicle  that  consists 
of  a  straight  track  and  a  maneuvering  track  and  is  shown  in 
Fig.  39. 

The  x-direction  vs  time  of  range,  velocity,  acceleration, 
bearing,  bearing  rate  and  bearing  acceleration,  were  plotted 
and  shown  in  Figs.  40-45. 

For  the  scenario  considered,  an  aircraft  ‘  •  =  10  KM  moves 

at  100  m/sec  (200  knots) ,  can  maneuver  at  a  maximum  accelera- 
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tion  of  4g  (A  =  39.24  m/sec  ),  and  has  a  probability  of 
maneuvering  at  max  0.2  (2P^  =  0.2,  P^  =  0.1),  and  a  probability 
0.5  of  not  maneuvering  at  all  (P£  =  0.5).  Assume  a  Lazy 
maneuver  that  will  provide  correlated  acceleration  inputs  for 
periods  between  10  and  30  sec.  Hence,  with  an  average  maneuver 
duration  of  20  sec,  by  using  Eq.  (4.19),  y  =  1/20  =  0.05  =  a, 
and  the  correlation  coefficient  P=l-aT=l-  .0G5T. 

The  two  sensors  were  classified  as: 
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FIG.  40.  RANGE  (X-AXIS)  VS  TIME 
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PIG.  41.  VELOCITY  (X-AXIS)  VS  TIME 
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FIG.  42.  ACCELERATION  (X-AXIS)  VS  TIME 
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FIG.  43.  BEARING  (X-AXIS)  VS  TIME 
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FIG.  44.  BEARING  RATE  (X-AXIS)  VS  TIME 
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FIG.  45 


BEARING  ACCELERATION  (X-AXIS)  VS  TIME 


A.  EXAMPLE  1  -  AIR  SEARCH  RADAR 


The  radar  data  rate  was  ten  samples  per  second  (T  =  10 
sec)  hence,  the  correlation  coefficient  p  =  1  -  0.05(10)  =  0.5, 
and  the  sensor  processing  noise  (measurement  noise  variances) 
has  been  taken  into  account,  500  m  in  range  (aD  =  500  m)  and 
17.4  mrad  in  bearing  (o.  =  1  degree  =  17.4  mrad) .  The 

U 

variances  of  maneuvering  at  max  acceleration,  and  not  maneu¬ 
vering  at  all,  were  calculated  by  using  Eqs.  (4.17-4.18) 

2  2  -4-2 

and  found  to  be,  =  46193  m/sec  and  a =  1.346  *  10  "  sec 

1 .  Kalman  Filter  Evaluation 

This  model  can  be  used  with  the  Kalman  filter 
Eqs.  (4.29-4.31-4.33)  to  determine  a  set  of  filter  gains. 

For  the  model  assumed,  by  using  Eqs.  (4.5,  4.11-4.14, 
4.25),  the  following  matrices  were  obtained: 
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2 •  Simplified  Kalman  Filter  Evaluation 

The  simplified  Kalman  filter  is  achieved  by  assuming 
(incorrectly)  that  the  vehicle's  change  in  velocity  is 
uncorrelated  between  samples;  i.e.,  the  maneuver  is  white,  and 
p  =  0.  Here,  for  the  model  assumed,  by  using  Eqs.  (4.2- 
4.4-4.5-4.26-4.27),  the  following  matrices  are  obtained. 
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EXAMPLE  2  - 

-  SURFACE 

AND  AIR 

SEARCH  RADAR 

In  this  example/  all  quantities  except  the  following/  have 
previously  been  defined  for  the  air  search  radar  example: 
Sampling  time  T  =  1  sec,  hence  p  =  1-aT  =  1-0.05(1)  =  0.95, 
the  sensor  processing  noise  (measurement  noise  variances) 
has  been  taken  into  account,  =  20  m,  oQ  -  0.1  degree  =  1.74 
mrad  and  the  variances  of  maneuvering  were  calculated  and 
found  to  be,  g3^  =  461.93  m/sec^,  g^2  =  1.346*10  ^  sec  3 . 

1 .  Kalman  Filter  Evaluation 

For  tne  model  assumed  the  following  matrices  were 
obtained* 
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C.  THE  ct-8  FILTER  EVALUATION 

The  o-8  filter  considered  in  this  paper,  and  used  for 
simulation,  is  one  of  many  varieties  possible  in  this  class, 
is  more  easily  implemented  than  either  the  Kalman  or  Simpli¬ 
fied  Kalman  filters,  and  has  been  selected  for  evaluation 
since  it  is  utilized  extensively  in  tactical  applications. 
Because  it  is  designed  to  minimize  the  mean  squared  error  in 
filtered  position  and  velocity  under  the  assumption  of  straight 
line  target  motion,  it  has  little  capability  to  track  severely 
maneuvering  vehicle.  The  a- 8  filter  examined  has  no  provision 
to  adapt  to  different  target  types,  as  does  the  Kalman  filter, 
since  maneuver  statistics  are  not  taken  into  account.  The 
Eqs.  (2.1-2.3-2.37-2.38)  for  the  a-8  filter  evaluated  are: 


3^  = 


N  N  N 
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i  X  -  4', 
4  -  + 
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where 

2 (2N-1 )  „  6  „  „ 

“  =  5Tn+1)~  '  8  =  STn+IT  and  T  =  1  sec  or  10  sec 

in  order  to  continually  compute  the  least-squares  line  through 
the  observations. 

D.  COMPARISON  OF  FILTER  ACCURACIES 

One  hundred  (100)  Monte  Carlo  trials  were  made  for  each 
combination  of  tracking  filter,  tracking  sensor  and  data  entry 
procedure. 

Experimental  determined  filtered  and  predicted  accuracies 
in  vehicle  range  coordinates  (range,  range  rate,  and  range 
acceleration)  and  bearing  coordinates  (bearing,  bearing  rate, 
and  bearing  acceleration)  with  means  and  variances  of  estima¬ 
tion  error  histories,  were  then  calculated  and  plotted  and 
shown  in  Appendix  A. 

Table  II  shows  a  representative  summary  of  the  result, 
in  which  the  prediction  accuracies  of  each  filter  were 
compared  on  a  percentage  basis  to  that  of  the  Kalman  filter. 

Tne  entries  in  the  table  were  determined  by  averaging  the 
experimentally  obtained  percentage  degradations  in  each  of 
the  above  coordinates. 

The  Simplified  Kalman  filter,  and  the  Kalman  filter 
generally  performed  within  twenty  percent  (20%)  of  each  other. 

The  a-B  filter  performance,  on.  the  average,  appears  to 
be  about  equal  to  the  above  filters  for  the  straight  part 
of  track,  and  about  thirty  to  sixty  percent  (30-60%)  worse. 


Key:  1  =  within  10-20  percent  of  the  Kalman  filter 
2  =  within  30-50  percent  of  the  Kalman  filter 

Table  II.  Synopsis  of  the  Accuracy  Comparison 
of  the  Three  Tracking  Filters 

with  the  greatest  degradation  occurring  for  the  maneuvering 
and  accelerating  part  of  track,  because  the  gain  vector  quickly 
becomes  too  small  to  correct  for  the  large  estimation  errors 
resulting  from  target  maneuvers. 


VI.  SUMMARY  AND  CONCLUSIONS 


While  no  extensive  analysis  is  implemented,  it  is  con¬ 
sidered  that  a  reasonable  and  unequivocal  comparison  of  the 
filters  can  be  made  from  the  material  presented. 

The  analysis  of  the  filters  and  the  radar  system  simula¬ 
tions  presented  in  this  paper  is  considered  overall  to  be  a 
simplified  but  realistic  model  of  a  sophisticated  system  which 
could  be  implemented  with  current  "state  of  the  art"  hardware. 
Based  on  the  ensemble  averages  the  Kalman  filter  obviously 
provides  a  somewhat  better  tracking  response  for  the  target 
track  tested. 

The  tracking  ability  of  the  a-g  and  Kalman  filters  appears 
to  be  about  equal  for  "look  alike"  targets  in  close  proximity, 
under  the  assumption  of  straight  line  motion.  The  a-g  filter, 
however,  provided  unsatisfactory  performance  when  the  tracked 
vehicles  executed  maneuver^ .  Based  also  on  the  simulation 
results,  the  Simplified  Kalman  filter  becomes  attractive  for 
implementation,  because  it  provided  tracking  accuracies  within 
ten-twenty  percent  of  the  Kalman  filter.  The  utility  of  this 
filter  is  greatest,  when  either  the  sensitivity  of  tracking 
performance  to  assumed  maneuver  correlation  is  small,  or  when 
the  target  maneuver  approaches  whiteness  relative  to  the 
sensor  data  rate. 

The  filter  implementation  requirements  increase  in  the 
following  order:  a-g  filter.  Simplified  Kalman  filter. 


Kalman  filter.  Moreover,  the  "complexity  factor"  between  the 
above  filters  is  about  two-to-one. 

Finally,  in  most  applications,  the  answer  to  the  question, 
"which  filter  is  most  accurate?",  does  not  alone  determine 
filter  selection.  Indeed,  the  following  questions  must  all 
be  answered  in  the  filter  selection  process  to  obtain  the 
"best"  filter  for  a  particular  system: 

a.  What  are  the  actual  accuracies  of  each  filter? 

b.  What  are  the  relative  filter  accuracies? 

c.  What  are  the  tracking  accuracy  requirements  of  the 
system? 

d.  How  sensitive  is  system  performance  to  tracking 
accuracy? 

e.  Whac  are  the  computer  requirements  of  the  filter? 

f.  What  are  the  computer  limitations  of  the  system? 

The  list  shows  thac  filter  selection  involves  careful  balancing 
of  filter  accuracies,  filter  implementation  requirements, 
and  system  performance  goals  and  limitations. 
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APPENDIX  A 


Simulation  results,  of  experimental  filtered  and  predicted 
accuracies,  m  vehicle  range  and  bearing  coordinates,  with 
means  and  variances  of  estimation  error  histories,  provided 
for  comparison  of  each  of  the  three  filters. 

These  results  were  generated  using  Monte  Carlo  simulation, 
with  the  following  description  of  run: 
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FIG.  46.  RANGE  VS  TIME  FOR  <x-6  FILTER 
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FIG.  47.  VELcJlTY  VS  TIME  FOR  o-8  FILTER 


Jit’ 


u.o 


7.  50 


15.  CC 


22.50 


30.00 


X-Scale:  0.375  units 
Y-Scale:  1.060  units 


FIG.  50.  VARIANCE  OF  RANGE  ERROR  VS  TIME  FOR  o-S  FILTER 
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FIG.  53.  VELOCITY  VS  TIME  FOR  SIMPLIFIED  KALMAN  FILTER 
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FIG.  54.  BEARING  VS  TIME  FOR  SIMPLIFIED  KALMAN  FILTER 
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FIG.  53.  BEARING  VS  TIME  FOR  SIMPLIFIED  KALMAN  FILTER 
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FIG.  57.  MEAN  VELOCITY  ERROR  FOR  SIMPLIFIED  KALMAN  FILTER 
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FIG.  59.  MEAN  BEARING  RATE  ERROR  VS  TIME  FOR  SIMPLIFIED  KALMAN 

FILTER 
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.  60.  VARIANCE  OF  RANGE  ERROR  VS  TIME  FOR  SIMPLIFIED 
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FIG.  63.  VARIANCE  OF  BEARING  RATE  ERROR  VS  TIME  FOR 

SIMPLIFIED  KALMAN  FILTER 
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FIG.  65.  VELOCITY  VS  TIME  FOR  KALMAN  FILTER 
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FIG.  66.  ACCELERATION  VS  TIME  FOR  KALMAN  FILTER 


(RAD/SEC) 


3.3001-02 


2.  2016-02 


*4  C25E-C3 


-J.  S43E-C3 


-l.«5*E-C2 


-2.«*E  a 


- - - - - - . » 


X-Scale:  0.375  units 
Y-Scale:  0,0013  units 


3.5C0E-O2 


2.  2C1H-02 


3.  0256-03 


-3.  '"*36-03 


-U«'3E-C2 


-2.*»3*E-02 


»  -*.2*36-02 

**,S0  30.00 

(SEC) 

( - ) - )  Actual 

( *  ♦  •  • )  Filtered 


FIG.  68.  BEARING  RATE  VS  TIME  FOR  KALMAN  FILTER 
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FIG.  69.  BEARING  ACCELERATION  VS  TIME  FOR  KALMAN  FILTER 
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FIG.  70.  MEAN  RANGE  ERROR  VS  TIME  FOR  KALMAN  FILTER 
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FIG.  71.  MEAN  VELOCITY  ERROR  VS  TIME  FOR  KALMAN  FILTER 
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FIG.  72.  MEAN  ACCELERATION  ERROR  VS  TIME  FOR  KALMAN  FILTER 
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PIG.  78.  VARIANCE  OF  ACCELERATION  ERROR  VS  TIME  FOR  KALMAN  FILTER 
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FIG.  80.  VARIANCE  OF  BEARING  RATE  ERROR  VS  TIME 

FOR  KALMAN  FILTER 
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